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C^l ' Under the conditions of weak Langmuir turbulence, a self-consistent wave-particle Hamiltonian 

models the effective nonlinear interaction of a spectrum of M waves with N resonant out-of- 
equilibrium tail electrons. In order to address its intrinsically nonlinear time-asymptotic behavior, 
a Monte Carlo code was built to estimate its equilibrium statistical mechanics in both the canonical 
^— I ' and microcanonical ensembles. First the single wave model is considered in the cold beam/plasma 

OQ , instability and in the O'Neil setting for nonlinear Landau damping. O'NeiPs threshold, that sep- 

arates nonzero time-asymptotic wave amplitude states from zero ones, is associated to a second 
order phase transition. These two studies provide both a testbed for the Monte Carlo canonical 
and microcanonical codes, with the comparison with exact canonical results, and an opportunity 
to propose quantitative results to longstanding issues in basic nonlinear plasma physics. Then the 
properly speaking weak turbulence framework is considered through the case of a large spectrum of 
waves. Focusing on the small coupling limit, as a benchmark for the statistical mechanics of weak 
Langmuir turbulence, it is shown that Monte Carlo microcanonical results fully agree with an exact 
microcanonical derivation. The wave spectrum is predicted to collapse towards small wavelengths 
Oh' together with the escape of initially resonant particles towards low bulk plasma thermal speeds. 

C/j , This study reveals the fundamental discrepancy between the long-time dynamics of single waves, 

O ' that can support finite amplitude steady states, and of wave spectra, that cannot. 

• i-H 
CZ2 . 

PACS numbers: 05.20.-y, 52.35.-g, 52.35.Ra, 05.10.Ln, 45.50.-j 

I. INTRODUCTION 

• Wave-particle interaction is a universal phenomenon in space and laboratory plasmas. Its importance ranges from 
magnetic confinement fusion devices, in particular with its relation to anomalous transport or supplementary heating, 
to laser plasma interaction, astrophysics and charged particle beam physics. It is responsible for high-frequency plasma 
turbulence which differs substantially from low-frequency fluid or magnetohydrodynamic turbulence. In this article, 
we shall consider one of the simplest wave-particle interaction settings, namely electrostatic wave-particle interaction 
■ taking place in bump-on-tail or weak beam-plasma instabilities. 

Under the conditions of weak Langmuir turbulence electrostatic wave-particle interaction may be modeled 
0, H, H, H, @ by a reduced system coupling self-consistently a set of M longitudinal (Langmuir) waves to N quasi- 
t^j ' resonant tail particles. This approach is especially appropriate to the description of basic plasma kinetic phenomena 
such as the bump-on-tail instability [7] or Landau damping |8|. In this reduced framework, the effective dynamics 
^sj' is ruled by a (M + ./V)-dimcnsional Hamiltonian system. More specifically, particles in the background plasma 
(whose velocities are roughly smaller than the thermal velocity) only participate to the dynamics via a small subset 
Oh 1 of collective modes, the Langmuir (long-wavelength) modes of longitudinal oscillation around their guiding-center 
' positions. These modes are described by action-angle variables (Ij,6j) with 1 < j < M for the M waves. In the 
absence of resonant particles, they oscillate with constant angular frequencies d9j/dt = u>oj, that, according to the 
Bohm-Gross dispersion relation, are approximately equal to the plasma frequency lj p for long-wavelength modes. 
These waves strongly interact with those plasma particles in the tails of the distribution function having velocities 
v close to ujQj/kj. The coupling is controlled by the small parameter rj that is the ratio of the tail density over the 
bulk plasma density. The dynamics of N identical quasi-resonant particles moving on the interval of length L with 
periodic boundary conditions, with unit mass and charge, and respectively position x r and momentum p r , interacting 
with M waves with wavenumbers kj = j2n/L, derives then from the Hamiltonian 
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Equations of motion are then 

xi = Pi, (2) 

f^T) M 32 

—L^cosikjXi-Oj), (4) 
3 i=i 

sin (kjxi — 9j) . (5) 

i=i 

It is important here to insist on the principle of this derivation: the background plasma is assumed to respond linearly 
to the waves which is valid provided resonant particles, whose velocities are much larger than plasma thermal velocity, 
only form a small fraction, t], of the total plasma. Moreover it is interesting to note that the same kind of wave-particle 
reduction has been undertaken for other physical regimes of wave-particle interaction. In particular, a Hamiltonian 
self-consistent wave-particle model Q has recently been built in order to study the nonlinear interaction of a packet 
of waves with a nonequilibrium electron distribution in a magnetized background plasma. 

Under the above hypotheses, the Hamiltonian model (Q} contains the effective dynamics of weak Langmuir tur- 
bulence. Yet, if a Hamiltonian system exhibits some ergodicity, then thermodynamics states that time averages 
can be replaced with equivalent space averages over the microcanonical ensemble. A Monte Carlo code was built 
to estimate the equilibrium statistical mechanics of (|TJ) in both the canonical and microcanonical ensembles. This 
provides an efficient tool to investigate the time-asymptotic fate of the wave-particle model and discuss the evolution 
of weak Langmuir turbulence. This is of particular interest in the broad spectrum case where direct thermodynamical 
derivations are difficult. 

In Section HH we present the canonical ensemble of the M waves/ N particles self-consistent Hamiltonian model and 
focus on the M = 1 case. In Sec. IIII1 we consider the cold beam/plasma instability where the single wave approxima- 
tion is particularly appropriate. Then, in Sec. IIV[ we apply the single wave model to the original O'Neil setting 
for nonlinear Landau damping and, in parallel with microcanonical Monte Carlo results, associate O'Neil's threshold 
to a second order phase transition. These two studies provide both a testbed for the Monte Carlo canonical and 
microcanonical codes, with the comparison with exact canonical results, and an opportunity to propose quantitative 
results to longstanding issues in basic nonlinear plasma physics. It is shown that statistical mechanics is a useful quan- 
titative tool to address the nonlinear time-asymptotic saturation stage. Quantitative comparisons with experiments 
are given. Then, in Sec. [Vj we turn to the properly speaking weak turbulence framework by considering the case of a 
large spectrum of waves. We focus on the small coupling limit, as a benchmark for the statistical mechanics of weak 
Langmuir turbulence. We show that Monte Carlo microcanonical results fully agree with an exact microcanonical 
derivation. We conclude by discussing the predicted time-asymptotic spectrum, its relation to quasilinear theory and 
to a pending debate in theoretical, experimental and numerical plasma physics on the long time fate of the electric 
field for Landau damping. 

II. CANONICAL ENSEMBLE 

We shall start deriving statistical mechanics in the canonical ensemble and then use the equivalence between 
canonical and microcanonical ensembles. Apart from the energy, it can be easily checked that the average total 
momentum per particle, a, is also conserved with 



a = P+jjJ2 k 3 I 3> ( 6 ) 

and 

1 N 

1=1 

We introduce the new particle momenta pi = pi — p and the reduced wave intensities 



*i = |. (8) 
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that are properly normalized with respect to the N — > oo limit [TTf| . Then, introducing ([6]), ([T]) reads 
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The canonical partition function is Z C (N,(3) = Zq(N, (3)Zi(N, pi, a), with 
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Z x (N, M, 13, a) = N M L N / J[ diP j d9 J exp [-N/3g(/3, a, ip, 9)] , 
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All the difficulty lies in the estimation of (|T3|) for M > 1 . We shall then first restrict our discussion to the case where 
only one mode is selected before going to the large-M case in the small-coupling limit in Sec. [V] The M — 1 case 
occurs naturally in the cold beam-plasma interaction 0, OH [H, [IH ■ We get, for k = 1 and L = 2ir, 



1 



2tt 



This gives 



709, a, ^ 0) = i (<r - V) 2 + woV> - 1 In [/ (/9^ /2 



(14) 
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Note here that in this particular case g{(3, a, ip, 9) is actually independent of 9. To perform the ip integral, we use 
Laplace's method, for which we need to find the value ip* of ip for which g(/3, a, ip, 9) becomes a minimum. We get 



9 9 

— = ip + uj q - a 
dip 



for 



ip* = a - uj + t]P3ujq 



(16) 



(17) 



We have, Vip > 0, dg/dtp > dg/dip\ — u n — <r — t]I3uIq/2. If ujq — a — t]/3uj^/2 > 0, g attains its minimum at the 
boundary, namely in ip = 0. Otherwise, there exists a unique minimum in ip* > satisfying (|17j) . Transition occurs 
for the critical 
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We have then 

/>00 

Z x {N,l3,cr) = 2irN / exp [-Nf3g(f3, a, f)] 
Jo 

/•OO 

~ 27riVexp[-iV/3 5 (/3,cr,^*)] / d^exp 

Jo 



The ensemble average of the wave intensity is then simply 

(I) c {N,p,v) = Nil>*(p,a)+o(N), (19) 

meaning that the wave intensity is extensive when tp* > 0, and non-extensive when ip* — 0. The mean intensity I/N 
acts as an order parameter (see Ref. Q): If f3 c < 0, then tp* is always positive; yet, if (3 C > 0, then for (3 < (3 C (high 
temperature regime), "0* = meaning that, as t — ► oo, the wave intensity ip should vanish and particle motions should 
be ballistic, whereas for (3 > [3 C (low temperature regime), ip* > meaning that the wave intensity "0 should be finite 
as t — > oo. 

III. THE SINGLE WAVE CASE: SATURATION OF THE COLD BEAM/PLASMA INTERACTION 

A. Presentation 

The analysis of O'Neil and Drummond et al. [l2l . [l3| have shown that, in the case of a monokinetic beam, only 
a single wave, the wave having the largest growth rate, has a significant interaction with the beam, at least up to 
the first trapping oscillations. We shall consider here this single wave case and retain one mode (M = 1) with rest 
frequency loq — ui p and wave number k = 1. A simple linear analysis [l2l . [l5| determines then which initial beam 
streaming velocity Vb leads to the maximal destabilization of this mode. This gives Vb = (l + 2 -4 / 3 7/ 1 / 3 ) u> p and 
associated maximal linear growth rate 7l = 2 _4 / 3 v / 3f? 1 ^ 3 Wp- This derivation is analogous to considering the more 
realistic scenario where the initial cold beam velocity is given and determining the dominant mode, namely the mode, 
among a continuum of modes, that is the most strongly destabilized. The expression of Vb reveals the weakly resonant 
nature of the linear instability: it concerns mostly particles having velocities not equal to the initial wave velocity but 
about 7l/& away from it [16j . Fig. [1] shows the time evolution of the wave intensity starting from an infinitesimal 
value. Nonlinear effects, namely trapping of particles in the wave trough, stop the initial linear exponential growth. 
Trapping oscillations coincide with the existence of a strong density inhomogeneity in the phase space. This is visible 
on the phase space snapshot in Fig. 
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B. Dynamical and Monte Carlo simulations; Canonical approach 

We shall use the M = 1 cold beam instability ,J% as a simple testbed for the Monte Carlo code. Let us then 
first derive the canonical ensemble average of the wave intensity in this case and look for an analytical estimate as a 
function of the problem parameters and initial data. 

The canonical ensemble average of energy is 

U = -d fi In Z c = -dp In Z (N, /?) - dp In Z\(N, /3, a), (20) 

with 

-dp\nZ {N,(3) = ^±, (21) 

and 



- dplnZ^P,*) = Ng(P,a,r) + o(N), 



(22) 
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FIG. 1: Wave intensity ip — I/N as a function of time (in lo^ 1 units) in the cold beam case for N = 10000 and r\ — 3.85 x 10~ 5 
(plain line). For the same parameters, ensemble averages of ip as a function of the number of Monte Carlo iterations for the 
microcanonical ensemble (bold line) and for the canonical ensemble with U — Eq (dots) (See Sec. MI Bl) . 
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FIG. 2: Snapshot of the one-particle phase space in the wave frame for t — SOOOtJjJ" 1 . Instantaneous separatrices of the wave 
resonance are shown. Dots represent the N = 10000 particles, which were initially distributed uniformly in x with the same 
velocity Vb- 



This gives, with e = U/N, 



(23) 



that has to be solved together with ((T7|) . As given above, we have a — lu p = crf^ujp with c = 2 4 / 3 . 
Let us look for solutions such that = y / 2i] r ip*ujp^ 2 = arj a uj p . Then (jTTJ) yields 



= cri 1/3 uj p 



al (Par) a u>^ 



(24) 



A maximal ordering in is obtained for f]ai] a ujp finite, ft gives f — a — 1/3 = 2a— 1, that is a = 2/3. It is important 
to note that this corresponds to the trapping scaling for the wave intensity saturation, namely the asymptotic trapping 
frequency, scales as the linear growth rate: LJb* ~ rf^uip ~ Jl- fntroducing the finite quantity d = f3a-q 2 ^ui 2 , we 
have then 
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Let us now assume that the canonical average of the energy density, e, is close to its initial value, £o — Eq/N. 
Because the wave intensity is initially vanishingly small, the identification e = So gives for the cold beam case under 
consideration 
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This gives finally 



1 - 21 ; J ° ( ^ ^ ) , 



to be solved together with (|23|) . We get 

a = 2 2 / 3 [Jl(d) + F 2 (d)} 1/3 = c- 1 (F 1 (d) + 2F 2 (d)) 



Using c = 2 4 / 3 , this yields 
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This gives numerically d ~ 1.56 and a ~ 1.313, yielding, with (pM)) and in the limit A~ — > oo, 



O.8677 



1/3,, 



•0(1). 



(31) 



Dynamical symplectic Q simulations indicate that, after some trapping periods, the average wave intensity I/N scales 
approximately as 0.78r; 1//3 [l^]. For the parameter 77 chosen in Fig. [TJ this gives (I/N) ~ 0.026 as AT — > 00 which is 
indeed in agreement with the time average of tp. The microcanonical Monte Carlo results plotted on the same Figure 
show an excellent agreement with this time average. Canonical Monte Carlo results assuming U = Ne$ is in fine 
agreement with the analytical result (|3ip . This is a practical way to bypass the microcanonical or time- asymptotic 
evaluation of (3 and relate directly canonical results to initial data. Here this result shows a ten-per-cent discrepancy 
with the exact microcanonical result. Moreover, we can equivalently express Eq. (|31[) in terms of the ratio between 
the asymptotic trapping pulsation and the linear growth rate yielding 
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1.67. 



(32) 



This is in agreement with the results of an early experiment by Mizuno and Tanaka [13] designed to test the single wave 
model. A monoenergetic beam is injected along a homogeneous magnetic field, making the system unidimensional, 
into a plasma. In this electron gun, the mechanism of the nonlinear saturation of a nearly single wave is reported 
through the first measurement of the trapping of beam electrons. The approximate experimental value, given in Ref. 
pjj, of the ratio Ub/lL associated to the first trapping pulsation is 0.83. As evidenced by Fig. [l] this can here be 
used as a fairly good approximation of the asymptotic value u>b*/7L- Yet, proceeding to a direct measure from the 
data plotted on the Figure 1 of Ref. [l7| . we obtain a somewhat larger value for this quantity, slightly above 1.2 and 
closer to the estimate (p?2"|) . 

Finally, it is interesting to note that an explicit low-dimensional modeling of the dynamics, along the lines of Ref. Q 
stressing the phase space division between the macroparticle and chaotic sea components, has been recently proposed 
by Antoniazzi et al. in the saturated regime of a single pass free electron laser around perfect tuning [1 SI ] - This 
model uses an Hamiltonian formulation that is a direct counterpart of the present Hamiltonian single-wave particle 
framework. 



IV. LANDAU DAMPING OF A SINGLE WAVE: PHASE TRANSITION AND O'NEIL'S THRESHOLD 



In Ref. [10|, O'Neil considered the damping of a single wave interacting effectively with a subset of resonant 
particles. He proposed a qualitative dynamical threshold conditioning the time asymptotic behavior of the wave: if 
the linear damping time |7l| _1 is far smaller than the initial nonlinear trapping time t&o of the particles in the wave 
trough, then the wave should not have enough time to experience nonlinear effects and should be completely damped 
as t — ► oo. On the contrary, if the nonlinear time t^q is smaller than |7l| _1 , the linear one, nonlinear effects (trapping) 
should forever control the dynamics. This can be already inferred from a linear analysis. As already apparent in the 
cold beam case, the particles controlling the linear regime are those quasi-resonant particles located at about 7l/& 
from LUp/k (l6| : if those particles are trapped in the initial wave (with half resonance width 2u!bo) right away, then the 
linear regime should be bypassed. However, this result is qualitative. Actually, in spite of major recent advances in 
the analytical and computational fields 0, HH El- 0, I24I l25l . [2(| , mostly in the full Vlasov-Poisson framework, 
as well as in experiments [27l |. no definite answer presently exists to this intrinsically nonlinear and time-asymptotic 
issue, in the sense that no quantitative threshold separating nonzero time-asymptotic wave amplitude states from zero 
ones is known. This is basically due to the facts that analytical reductions of the problem should e.g. make some 
restrictions on the asymptotic behavior of the electric field whereas numerical and experimental investigations should 
be obviously taken with great care when addressing the t — > oo limit. For instance, the nice experiments recently 
done in a trapped pure electron plasma [27| diagnose the first trapping oscillation frequency uif, at the onset of the 
nonlinear stage (and more precisely the dimcnsionlcss ratio jL/uJb) as a function of the initial ratio jL/wbo rather 
than the, rigorously speaking, time- asymptotic state. 

We shall consider here the single wave-particle interaction model that describes the effective (field envelope) dynam- 
ics of the Vlasov-Poisson system in the case where only one wave exists. This is obviously an idealized or academic 
situation, yet it is perfectly appropriate to discuss O'Neil's picture. We shall show here that O'Neil's threshold can 
be related to a second order phase transition for the single wave self-consistent Hamiltonian model H ' , with order 
parameter the wave intensity if). This was already discussed in Ref. @. We shall however consider here specifically 
the O'Neil setting, that is we take a given finite initial wave intensity i/'o an d a monotonously decreasing, linear, tail 
particle distribution function fo(v), of which we shall vary the slope and thus the linear Landau damping rate 7l 
(See Fig. [3]). Because of locality properties in the velocity space |3Q|, the wave has an effective interaction only with 
particles that are not too far away from its separatrices. One can reasonably think that particles having velocities 
more than, let's say, icobo/k away from uj p /k [30] should evolve on Kolmogorov-Arnold-Moser (KAM) tori for the 
pendulum-like 1,5 degrees of freedom Hamiltonian (given by Eqs. ©-© with (I(t),9(t) j) associated to H N1 , mate- 
rializing the separation between the wave-particle interaction zone and the plasma bulk captured into the collective 
variable /. 
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FIG. 3: Family of tail distribution functions. We shall take A = Aubo/k. 

Let us then put ui p = 1, define A = iuibo = 4 (2rjipo) 1 ^ i and take fo(v) = av + b if v 6 [1 — A; 1 + A] and 
fo(v) = otherwise, with a < 0. Normalization imposes b = (1 — 2aA)/(2A) and positivity —a < (2A 2 ) -1 . We have 
(v) Q = 2aA 3 /3 + 1 and T = (t; 2 ) - (vf = -A 2 (4a 2 A 4 - 3) /9. Let us now write the condition (3 < (3 C for which 
particle asymptotic behavior should be ballistic. From ([23]) . we have in this regime 
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If we identify, once again, this canonical ensemble average with the initial density of energy in the system Eq — 
(v 2 ) /2 + ip and use a = (v) + t/jq, we get the following condition on the slope of the initial tail distribution: 

27 

f3 < (5 C Q(a) > - V A- 9 (34) 
16 

where Q is the third degree polynomial 

Q(x) = (x — c + ) (x — c_) (x — c) (35) 
with c± = c± V3/(2A 2 ),c= -3A/(4 5 ?7). (36) 

Moreover, to ensure (3 C > 0, one should have a < ujq i.e. (v) + ipo < 1 so that 

a<c= -3A/(4 5 r?). (37) 

Q'(x) vanishes for d± = c±l/ (2A 2 ), the relative maximum of Q being reached in d_. Taking into account the positivity 
condition —a < (2A 2 ) -1 , Eq. (|37p and d_ < — (2A 2 ) -1 , one concludes that, in order to observe a phase transition for 
this family of initial tail distribution functions, it is necessary and sufficient that Q(— (2A 2 ) -1 ) > 2777A -9 /16. The 
necessary condition Q > H^A -9 is easier to formulate and reads 

A > 3 x 2- 2 /y/ 3 . (38) 

The righthand side of (|38[) is reminiscent of the scaling law of the growth rate of the cold beam. This inequality is 
both a condition on the beam width and initial wave amplitude. In order to observe a regime in which the wave and 
the particles are decoupled in the sense that / is no longer extensive (ballistic regime instead of trapping regime) , it is 
necessary to exclude the case of a cold beam. In this latter case, even if the initial wave amplitude is infinitesimal, it 
will not be completely damped: this signals the singular character of the cold beam/wave interaction. Eq. (|38p means 
that a wave interacting resonantly with a cold beam cannot be Landau damped to zero. Additionally, the condition 
(|38|) thus specifies what is meant by an "initially finite" wave amplitude. The wave-particle interaction model a la 
O'Neil presented here rules out the classical research case of Landau damping in the limit of a perturbative wave 
of infinitesimal amplitude. The present framework is definitely non trivial since, if the wave is effectively damped 
asymptotically to zero, nonlinear effects should place the study in the frame of Isichenko's resolution |2(| predicting 
a very slow, asymptotically algebraic, damping. 

Let us illustrate these results for values of A and r/ allowing a phase transition in the canonical approach just 
described. Taking A = 0.25 and r) = 0.001 gives Q(-(2A 2 )- 1 ) = 1011.5 > %V&~ 9 = 442.37, which ensures 
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FIG. 4: Monte Carlo microcanonical ensemble averages of the wave intensity ip = I/N as a function of the slope a of the initial 
distribution function and corresponding q. For TV = 10000 particles, the microcanonical ensemble averages of the fluctuations 
< (ip— < if) > M ) 2 >i are indicated as errorbars. 

the existence of a phase transition. The tail distribution slope a can vary between —8 and 0, and the critical 
slope associated with (|34[) is a c = —3.106. This corresponds to a linear Landau damping rate jlc = (fta-cV) /2 = 
—4.88 x 10 -3 . We are now able to estimate the O'Neil threshold, namely the critical threshold of q = \ /^bo, as 
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0.078. 



(39) 
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FIG. 5: Same as Fig. [4] in lin-log scale. 



We can check that this value is very close to the threshold estimated in Ref. [8j in which the initial tail distribution 
(close to, but not, linear), and thus the linear Landau damping rate, was given whereas the initial wave amplitude 
was the control parameter. This agreement may be interpreted as an indicator of the robustness of the statistical 
mechanics approach. Figs. [4] and \5\ show the Monte Carlo microcanonical ensemble averages of the wave intensity 
as a function of the slope of the initial tail distribution function and of the parameter q. A nice agreement with the 
threshold ([59")) is found. The transition is emphasized by the logarithmic scale in Fig. as in the asymptotic ballistic 
regime, ifj scales as 1/N. 

Using equilibrium statistical mechanics allows then to set in a quantitative way O'Neil's threshold. However strong 
metastability effects and long relaxation times towards equilibrium in the vicinity of the phase transition [8] temper 
the practical interest of this result, especially in the purely collisionless Vlasovian N — ► oo limit [ll| for which there is 
increasing evidence [7l. l28l . l29l ] that it may not commute with the t — > oo limit. This explains that practical experimental 
and numerical attempts to determine q c typically yield values that are about one order larger than (|39|) . It should be 
also remembered that, as demonstrated by Isichenko's calculation [20| . the asymptotic behavior of a damping electric 
field is intrinsically very slow, so very difficult to detect and discriminate practically from a Bernstein-Greene-Kruskal 
(BGK) saturation stage [lij . Equilibrium statistical mechanics should prove however essential to address the large 
time nonlinear fate of wave-particle interaction in the generic case of many waves. Strong resonances overlap may 
then create large scale chaos inducing an efficient sweeping of the phase space ensuring good ergodic properties or, 
using the Lynden-Bell terminology, violent relaxation. 



V. MICROCANONICAL PREDICTIONS FOR A LARGE NUMBER OF WAVES: THE 

INFINITESIMAL COUPLING LIMIT 



A. Motivations and Monte Carlo simulations 

We shall consider here a regime for which the validity of quasilinear theory has been questioned [3ll . l32l . |33| , when 
the nonlinear timescale (resonance broadening time trb) is much lower than the linear one (of the order of the inverse 
of the linear growth rate 71,), that is fj, = (trb7l) 1 3> 1. Equilibrium predictions are expected to be particularly 
relevant for this problem since chaotic mixing and system ergodization should prove effective typically on the nonlinear 
timescale. 

Let us begin by presenting microcanonical Monte Carlo simulations done for parameters very close to those used 
by Doxas and Cary in Ref. [32]. Specifically, the wave spectrum is composed of M = 584 waves, that are initially 
strongly overlapping with a Chirikov overlap parameter of about 500, with wavenumbers between 2.78 and 3.70. The 
strong overlap is needed to model a continuous spectrum. The beam distribution in velocity is the same as in Ref. 
[IH, chosen to give a constant linear growth rate across the spectrum: fo(v) = C\ + C2 (1 — v a /v) if v G [w Q ;wb] 
and fo(v) = otherwise, with C\ and C2 fixed by normalization and the value of the growth rate. For the Monte 
Carlo code, these particular initial conditions just fix the two microcanonical constants of motion, namely the total 
energy and momentum. Then, as detailed in the analysis of Ref. [32j, the large- /1 regime is entered as the coupling 
(control) parameter 77 becomes vanishingly small. We are thus investigating here the small coupling regime. Results 
may be summarized as follows: in the transient stage of the Monte Carlo runs, that may be compared with some 
care to the real time dynamics, one observes "noisy" wave spectra that may be related to the spontaneous spectra 
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FIG. 6: (Color Online) Wave spectrum obtained through microcanonical Monte Carlo simulations, for r) — 2.6 x 10~ 4 and 
initial conditions described in the text, for two values of the number N of particles: with N = 5000 in the transient regime 
(after only a few iterations) (upper curve), with N = 5000 at equilibrium (middle curve) and for N = 20000 at equilibrium 
(lower curve). Wave amplitudes are given in logarithmic scale. 



discretization reported in previous studies [32]. In the thermodynamical framework, this appears only as a transient 
stage towards equilibration. The major point depicted by Fig. [5] is that the energy apparently tends to concentrate 
at equilibrium into a very small number of small wavelength modes while the long wavelength wave intensities are no 
longer extensive. 



B. Exact microcanonical derivation 



Let us consider analytically the small coupling limit, namely rj — > + . This is a singular limit that amounts to 
neglect (at zero order in rf) the wave-particle interaction potential in the self-consistent Hamiltonian ([I]) while retaining 
the total wave-particle momentum conservation constraint. In the microcanonical ensemble, the partition function 
reads, with uiqj = w p for any wave number j, 



. I N 2 M \ / N M \ 

\ 1=1 3=1 / \ 1=1 3=1 J 

Let us introduce Jj = Ij for 1 < j < M — 1 and I' M such that Y^=i Ij = Z^jLi 1 Ij + = MI' M . We get 

n / N 2 \ / N M ~l \ 

Z^E, P) = M J 5\E~Y^- uj p MI' M j s[P-^2m-J2 ^ ~ kM ^ r i ~ k MMI' M \ d N pd M I'. 



Introducing 



A/ 



. Im-i n \ 



(40) 



(41) 



B^E+pL.p, (42) 
w p 2 k M 

integration over I' M gives 

r I N M_1 i, u N z\ M_1 
Z,(E, P) = M MP-^-^fe- M 3_^25+^£& d"pJ] d /j 

J V J=l 3=1 ^ Wp 1=1 1 J = 1 

-m/, <*„ - W jj - b + £ - ^g) ] rfW-r 
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It is requested that 



M-l 

0< {k M -kj)I^<B, 



(44) 



so that, in particular, with ui p /kM = vg. 



< e. 



(45) 



This inequality (|45p is basically satisfied when the particle impulsions are larger than the smallest wave velocity, 
which is the natural setting of the wave-particle interaction model. Eq. (|43p is easily integrated over the pi 's yielding, 
assuming N even, 



Z^E,P) 



2M ( 2ituj v 



tv (N/2 - 1)! V k 



A I 



\ N/2 f I M-l 

) /#* 



(46) 



with f(x) = (B — x)^ N and where the integration boundaries are given < J^^Li 1 Q^M — kj) Ij < B with < Ij, 
Vj, 1 < j < M — 1. We begin by making a change a variable, putting (&m — kj) Ij = Ij. This gives 



2M 



tv (N/2 - 1)! V k 



2ivuj r 



< : m 



N/2 M-l 



'M-l 



Y[(k M - k j)- l [ / E^K 1 



(47) 



3=1 



We have 



M-l 



f\Yh\ d M ^i 



Jdlu-iJ... j dJi(s-^I 



(N-l)/2 



3 = 1 

M-l j 



M-l 



(N-l)/ 

M-2 



(JV-l)/2+l 



n 



3 

M-l 

n 



l\ (N-l)/2 + jJ 
1 



LA (jV-l)/2 + j 



dlM-i (B — Im-i 

B (N-l)/2+M-l 



3=2 

(iV-l)/2+M-2 



Finally, we get 



2Af 



„ v AT/2 M-l , s-1 



7T (N/2 - 1)! V. fc M J f \ (N - l)/2 + j 



—E-\ p - -P 

Up 2k m 



-i (iV-l)/2+M-l 



(48) 



Following the calculations presented in the Appendix El we derive the wave intensities for any jo strictly smaller than 

M as 



/ / \W 

W3o) u 



i r M v ejo 



(49) 



N UJ p V 6 j - v 0min 

As shown on Fig. this result nicely agrees with Monte Carlo microcanonical expectations in the limit rj — > + . We 
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wave amplitudes 




FIG. 7: (Color Online) Equilibrium wave spectrum as derived in Eq. (|49p for the same initial wave and particle distributions 
as in Fig. [6] for two values of the number N of particles: for N = 5000 (upper curve) and for N — 20000 (lower curve). 
Wave amplitudes are given in logarithmic scale. Numerical data for the initial conditions are t>e m in = 0.2379, e = 0.0555 and 
a = 0.333. 



have introduced the microcanonical temperature T M through [3£ 



1 d In Z u 1 



T u dE 



2 km.. 



(50) 



2k * 



so that 



e + 



av 6l 



(51) 



This clarifies the inequality (|45[) . that just means that the temperature 7), should be positive. Let us now complete 
this study by moving to a grand canonical ensemble, that is both canonical with respect to the energy and the total 
momentum. 



C. Grand canonical approach 



We define 



Z G (/3, A) = J exp (-(3H - XP) d N pd M I, 



(52) 



nth H = Y^iLi i>~ + w p Sjii Ij an d P — YliLi Pi + Sfci kjlj. This gives immediately 



„ \ N/2 / at\2\ M 

Z r ,r>.\) ( jj exp(— ) Hi H • A/, v 1 



(53) 



In this ensemble, we easily derive 



with (see Appendix |B|) 



1/(3 



30/ G 



(ui p + k J0 \/(3) io p (V0j o + \/(3) ' 



(54) 



/3 \N 



1 N 
Iv 51 



Pi 



1=1 I G 



(55) 
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and the ensemble average of the reduced kinetic energy is found to be 

/„ i / i „ \ " , 



i G 



eM^-^x>) 

W=l \ k=\ J 



G 

\ 2, 



m,-m-) ).-r 



D. Discussion and conclusion 



Considering a spectrum of waves, we have focused on the small coupling limit as a benchmark for the statistical 
mechanics of weak Langmuir turbulence. A trivial correspondence, assuming the equivalence of ensembles, between 
Eqs. P5|) and using ([55]) . shows that electrons should, as time proceeds, escape from the original resonances 

towards lower speeds. Meanwhile, the wave spectrum collapses towards short-wavelengths and the wave energy 
eventually concentrates into the mode of minimal phase velocity, that basically matches plasma thermal velocity. 
Consistently, particles (electrons) behaving almost freely loose their resonant character and follow the decay of waves 
towards small velocities. 

Let us then briefly comment on the relation of these results with the investigation on the limits of quasilinear 
theory undertaken by Doxas and Cary [32J. It is obviously uneasy to draw conclusions on the out-of-equilibrium 
behavior, such as diffusion, from equilibrium statistical mechanics results. Some statements are still possible: When 
the nonlinear timescale is far larger than the linear one, quasilinear theory is expected not to apply as the timescale 
of validity of (almost) zero-order effects, such as the reaction of the resonant wave spectrum on the zero-order 
distribution function, becomes negligible. Actually, we infer from the equilibrium microcanonical results that, within 
some nonlinear thermalization timescale, the initial particle distribution function is substantially modified whereas 
the wave spectrum collapses towards the short wavelengths. The thermalization stage is expected to exhibit a noisy 
wave spectrum, as in Monte Carlo transients, as the wave spectrum collapse proceeds through higher order wave- wave 
couplings, which may explain the dynamical findings of diffusion enhancements [32| . 

Now, we can consider the small coupling results from a more general perspective. Let us consider an initial very weak 
beam-plasma or bump-on-tail instability destabilizing some spectrum of waves towards suprathermal levels in a wave- 
particle interaction dynamics captured by the self-consistent Hamiltonian |Tj). Then, the previous calculation predicts 
that this out-of-equilibrium dynamics cannot sustain itself as the wave energy should drop and cascade towards 
short-wavelengths with phase velocities reaching the plasma bulk thermal velocity. Obviously, the self-consistent 
Hamiltonian model reaches at this point its validity limit. This behavior is an effect of considering a spectrum of 
modes rather than single modes, as in Sees. IIIII and lIVl for which trapping ensures the long-standing stabilization of 
the instability. 

Recently a pending debate emerged concerning the time-asymptotic state of a large amplitude Landau damped 
wave |27j |. It is commonly expected that the system will eventually enter some stationary BGK equilibrium with a 
time-asymptotic finite wave amplitude provided the initial amplitude is large enough. This agrees with the statistical 
mechanics predictions for a single wave derived in Sec. IIVI under the O'Neil setting. The existence of BGK steady 
states is supported by some theoretical analysis and numerical results [2ll . I22I [23j . However, Brodin's numerical 
simulations [19| suggested that the wave amplitude never settles to a steady value as the energy of the resonant 
particles diffuses slowly into higher harmonics. Moreover Isichenko [20] argued that nonlinearities should not stop the 
damping of the electric field towards a vanishing time-asymptotic amplitude. Our point of view here is that, within 
the weak coupling hypothesis - that besides fully justifies the linear response of the background plasma -, the Landau 
damping of a spectrum of waves is predicted not to reach a BGK finite amplitude steady state with the resonant tail 
electrons drifting with wave energy towards plasma bulk. 



APPENDIX A: EXACT MICROCANONICAL ENSEMBLE AVERAGES OF THE WAVE INTENSITIES 

IN THE LIMIT 77 -> 0+ 



We wish to calculate 



N/2 r I M-l 




We have 



dZ, 



(Af+2) 



dk 0o n(N/2 + 



' ni (^A n/2+1 n + i f T , I M ^ (h , u , 



(JV-l)/2 



so that 



, lW fN t \N k M dzi N+2) {E,P) 



1 



2 / 2 (Ar + 1) dk j0 
and the microcanonical ensemble average of the jg wave intensity is 



N/2 8Zl N+2) /dk 10 



,< 1 ^27TWp^ ( JV , + 1 ) 



(iV) 



We have, from Eq. (|48j). 



2M / 2ttcj 



tt(JV/2 + 1)! V fcM 



W/2+1 



M-l 



(fc M - kj)~ 



— — ^- X 

1 (N + l)/2 + j 



E 



u> n 2k 



P 



M 



(N+l)/2+M-l 



(k M -kj ) 1 



so that 



, ){N ) = 2N (N-l)/2 + l 
^ J0 V (iV + 1) (N- l)/2 + M 

This gives finally, for jo ^ M, 



-J- n, 



k M LO p 

uj p 2k M 



x (k M - k jo y 



v 8jo ~ v 9min 



with T M the microcanonical temperature given in (|5f p . 



APPENDIX B: GRAND CANONICAL APPROACH: EXPLICIT RESULTS 



Let us consider 



Z G (/?, A, ax, a 2 , a 3 , a 4 ) = / exp I -/3ai y - /3a 2 w p Jj - Aa 3 - Aa 4 ^ kjlj J d^pa^I 

\ 1=1 3 = 1 1 = 1 3 = 1 j 

-, N M 



We have then immediately 



exp ^— /? a i~^ — Aa3pJ dp y exp (~f3a,2UJ p Ij — Xa&kjlj) dlj 



12m 



/3a i 



^ [N\ 2 a 2 \ M 

I exp ^ 2/3Q! 3 ) n (@ a2UJ P + ^kj) 1 



1 dlnZ G 



/ N 



11,2,3,4— 1 



iVA 

' /3 ' 
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and 



This gives 



f N 

E 



pi 



, d 2 Z G 
1 dai 



A 2 Z„ 



Ctl,2,3,4 



= 1 



'^P±\ l d\nZ G 



2 / 2AT 
(=1 / G 



Pi 



Y 2 A 2 TV 



/3 2 /?' 



_ A TV A 2 

11,2,3,4 — J- ^ r r~ 



A 1 TV 
^~2^~^' 



(B2) 



(B3) 



(B4) 
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